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Abstract 

A partial wave analysis of the centrally produced K + K~ and KgKg channels has been 
performed in pp collisions using an incident beam momentum of 450 GeV/c. An unambiguous 
physical solution has been found in each channel. The striking feature is the observation of 
peaks in the S-wave corresponding to the /o(1500) and /j(1710) with J = 0. The D-wave shows 
evidence for the /2(1270)/ei2(1320), the /2(1525) and the /2(2150) but there is no evidence for 
a statistically significant contribution in the D-wave in the 1.7 GeV mass region. 
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One of the fundamental predictions of QCD is the existence of glueballs. Current theoretical 
predictions based on lattice gauge calculations indicate that the lowest lying scalar glueball 
should be in the mass range 1500-1700 MeV 0. The / (1500) and the /j(1710) display clear 
glueball characteristics in that they are both produced in glue-rich production mechanisms and 
are either not seen, or are heavily suppressed in normal hadronic interactions. In addition, 
the /j(1710) is observed to decay dominantly to KK and yet it is not produced in K~p 
interactions M. However, the spin of the /j(1710) is still uncertain with J = or 2 being 
possible. 

In 1987 the MARK III collaboration published the results of a spin analysis of the /j(1710) 
observed in radiative J/ip decays 0. The analysis assumed that the /j(1710) region was pure 
spin zero or pure spin two and ignored interference effects. The results from the spin analysis 
showed that spin two was preferred over spin zero. 

In 1989 the WA76 collaboration published results from the analysis of the centrally produced 
KK system [|J. In an attempt to assess the spin of the /j(1710) the angular distributions in 
the 1.5 and 1.7 GeV mass regions were studied. The angular distributions in the two mass 
regions were found to be similar, and since it was assumed that the signal at 1.5 GeV was due 
to the /2(1525), it was concluded that the signal at 1.7 GeV was also spin two. From these two 
observations the /j(1710) was allocated spin two. 

A more sophisticated analysis of the Mark III data including both spin zero and two 
amplitudes and the possibility of interference between them showed that the 1.7 GeV mass 
region was dominated by spin zero. However, these analyses have only ever been published in 
conference proceedings and have never been followed up with publications in refereed journals. 

More recently, the E690 experiment at Fermilab has published the results of a partial wave 
analysis of the centrally produced KgKg system 0. The mass spectrum looks very similar to 
that of the WA76 experiment; however, the peak at 1.5 GeV is found to have spin zero and not 
spin two as had been assumed in the WA76 analysis. Unfortunately, above 1.58 GeV, due to 
an ambiguity in the solutions, a unique determination of the spin of the /j(1710) could not be 
made by E690. 

It is essential to know the spin of the /j(1710) because measurements of its production rate 
in J I if) decays indicates that if it has spin two then it would be consistent with being a normal 
qq meson whereas if it has spin zero it would be consistent with having a significant glueball 
component [0. 

This paper presents a study of the K + K~ and KgKg final states formed in the reaction 

VP -> Pf(KK)p s (1) 

at 450 GeV/c. The subscripts / and s indicate the fastest and slowest particles in the lab- 
oratory respectively. The WA102 experiment has been performed using the CERN Omega 
Spectrometer, the layout of which is described in ref. ||. 

The reaction 

pp-^p f (K + K-)p s (2) 

has been isolated from the sample of events having four outgoing charged tracks, by first 
imposing the following cuts on the components of the missing momentum: |missing P x \ < 14.0 
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GeV/c, |missing P y \ < 0.16 GeV/c and |missing P z \ < 0.08 GeV/c, where the x axis is along the 
beam direction. A correlation between pulse-height and momentum obtained from a system of 
scintillation counters was used to ensure that the slow particle was a proton. 

In order to select the K + K~ system, information from the Cerenkov counter was used. One 
centrally produced charged particle was required to be identified as an ambiguous K/p by the 
Cerenkov counter and the other particle was required to be consistent with being a kaon. The 
largest contamination to the real K + K~ final state comes from the reaction pp — > Aj + p s 7r~ 
where the high momentum tt + from the decay of the A ++ is misidentified as a K by the 
Cerenkov system. In order to reject this contamination the positive particle was assigned the 
7r mass and the A ++ (1232) signal has been removed by requiring M(pfir + ) > 1.5 GeV. 

The method of Ehrlich et al. ||, has been used to compute the mass squared of the two 
centrally produced particles assuming them to have equal mass. The resulting distribution is 
shown in fig. p]a) where peaks can be seen at the tt, K and p masses squared. A cut on the 
Ehrlich mass squared of 0.14 < M\ < 0.55 GeV 2 has been used to select a sample of 30 868 
K + K~ events. 

The peak at the 7r mass squared results from 7rs being misidentified as Ks by an inefficiency 
of the Cerenkov counter. To determine the contamination from these ir + n~ events inside the 
Ehrlich mass cut, real 7r + 7r~ events have been passed through a simulation of the apparatus 
in which the efficiency of the Cerenkov counter (94 ± 1 %) has been taken into account. The 
resulting distribution is shown as the shaded histogram in fig. |I]a). The amount of contamination 
is small and the effect this contamination has on the angular distribution has been found to be 



negligible |10 



The centrally produced K + K~ effective mass distribution is shown in fig. |T|b). The main 
features of the spectrum are evidence for a sharp threshold enhancement and peaks in the 1.5 
and 1.7 GeV regions. 

A Partial Wave Analysis (PWA) of the centrally produced K + K~ system has been per- 
formed assuming the K + K~ system is produced by the collision of two particles (referred to 
as exchanged particles) emitted by the scattered protons. The z axis is defined by the mo- 
mentum vector of the exchanged particle with the greatest four-momentum transferred in the 
K + K~ centre of mass. The y axis is defined by the cross product of the two exchanged particles 
in the pp centre of mass. The two variables needed to specify the decay process were taken as 
the polar and azimuthal angles (9, <ft) of the K~ in the K + K~ centre of mass relative to the 
coordinate system described above. 

The acceptance corrected moments, defined by 

I{Q) = Y.tmYl{n) + 2 £ t LM Re{Y^(Q)} (3) 

L L,M>0 

have been rescaled to the total number of observed events and are shown in fig. |2|. As can be 
seen the moments with M > 2 and L > 4 are small (i.e. £43, £44, and t^o) and hence only 
S, P, and D waves with m < 1 have been included in the PWA. The t Q0 moment represents 
the total acceptance corrected mass spectrum. As can be seen, by comparing this spectrum 
with the raw mass spectrum shown in fig. |I|b), the effect of the acceptance is to increase the 
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proportion of events with masses less than 1.3 GeV. This is mainly due to the effects of the 
momentum threshold of the Cerenkov counter. 



The amplitudes used for the PWA are defined in the reflectivity basis |frTJ . In this basis the 
angular distribution is given by a sum of two non-interfering terms corresponding to negative 
and positive values of reflectivity. The waves used were of the form with J = S, P and D, 
m = 0, 1 and reflectivity e = ±1. The expressions relating the moments (£lm) and the waves 
( J^J are given in table [TJ. Since the overall phase for each reflectivity is indeterminate, one wave 
in each reflectivity can be set to be real (Sq and for example) and hence two phases can be 
set to zero {4> s - and <p p + have been chosen). This results in 12 parameters to be determined 
from the fit to the angular distributions. 

The PWA has been performed independently in 40 MeV intervals of the K + K~ mass spec- 
trum. In each mass bin an event-by-event maximum likelihood method has been used. The 
function 

N 

F = ~ E ln{I(ty} + £ t LM e LM (4) 

i=l L,M 

has been minimised, where N is the number of events in a given mass bin, 6lm are the effi- 
ciency corrections calculated in the centre of the bin and ttM are t ne moments of the angular 
distribution. The moments calculated from the partial amplitudes are shown superimposed on 
the experimental moments in fig |2|. As can be seen the results of the fit reproduce well the 
experimental moments. 

The equations that express the moments via the partial wave amplitudes form a non-linear 
system that leads to inherent ambiguities. For a system with S, P and D waves there are 
eight solutions for each mass bin. In each mass bin one of these solutions is found from the 
fit to the experimental angular distributions; the other seven can then be calculated by the 



method described in ref. [11[]. In order to link the solutions in adjacent mass bins, the real and 
imaginary parts of the Barrelet function roots are required to be step-wise continuous and have 
finite derivatives as a function of mass |12| . By definition all the solutions give identical moments 



and identical values of the likelihood. The only way to differentiate between the solutions, if 
different, is to apply some external physical test, such as requiring that at threshold the S-wave 
is the dominant wave. 

The four complex roots, Zj, after the linking procedure are shown in fig. [TJc) and d). As 
can be seen the imaginary parts give little help in the linking procedure. However, the real 
parts are well separated in most places and hence it is possible to identify unambiguously all 
the PWA solutions in the whole mass range. In the 1.8 GeV mass region two roots do become 
close together, but if these two roots are swapped it results in events from the S-wave being 
transferred almost entirely into the P-wave and produces mass spectra that look unphysical. 
In addition, the zeros do not cross the real axis and hence there is no problem with bifurcation 
of the solutions. Near threshold the P-wave is the dominant contribution for five out of eight 
solutions, another is dominated by D-wave and another has the same amount of S-wave and 
P-wave. These seven solutions [|TUJ have been ruled out because the K + K~ cross section near 
threshold has been assumed to be dominated by S-wave. The remaining solution is shown in 
fig. §. 

The S-wave shows a threshold enhancement; the peaks at 1.5 GeV and 1.7 GeV are inter- 
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preted as being due to the / (1500) and fj (17 10) with J = 0. The D-wave shows peaks in the 
1.3 and 1.5 GeV regions, presumably due to the / 2 (1270)/a 2 (1320) and /2(1525) and a wide 
structure above 2 GeV. There is no evidence for any significant structure in the D-wave in the 
region of the /j(1710). In addition, there are no statistically significant structures in any of 
the other waves. 

Extensive Monte Carlo simulations have been performed to check the validity of this result. 
They show that the feed through from the S-wave to the D-wave is approximately 5 % from 
threshold to 1.3 GeV decreasing to less than 1 % for masses greater than 1.5 GeV. The feed 
through from D-wave to S-wave is found to be negligible over the entire mass range and hence 
gives confidence in the above observation that the /j(1710) has J = 0. 

A fit has been performed to the Sq wave using three interfering Breit-Wigners to describe 
the /q(980), /o(1500) and /j(1710) and a background of the form a(m — m t h) b exp(—cm — dm 2 ) , 
where m is the K + K~ mass, m t h is the K + K~ threshold mass and a, b, c, d are fit parameters. 
The Breit-Wigners have been convoluted with a Gaussian to account for the experimental mass 
resolution (a = 6 MeV at threshold rising to 19 MeV at 2 GeV). The resulting fit is shown in 
fig. Ha) and gives 



/o(980) M = 985 ± 10 MeV, r = 65 ±20 MeV 
/o(1500) M = 1497 ± 10 MeV, T = 104 ± 25 MeV 
/o(1710) M = 1730 ± 15 MeV, T = 100 ± 25 MeV 

parameters which are consistent with the PDG [|l]| values for these resonances. 

A fit has been performed to the Dq wave above 1.2 GeV using three incoherent relativistic 
spin 2 Breit-Wigners to describe the / 2 (1270)/a 2 (1320), / 2 (1525) and the peak at 2.2 GeV a 
background of the form a(m — m t h) b exp(—cm — dm 2 ), where m is the K + K~ mass, m t h is the 
K + K~ threshold mass and a, b, c, d are fit parameters. The resulting fit is shown in fig. |]b) 
and gives 

/ 2 (1270)/a 2 (1320) M = 1305 ± 20 MeV, T = 132 ± 25 MeV 
/£(1525) M = 1515 ± 15 MeV, T = 70 ±25 MeV 

/ 2 (2150) M = 2130 ± 35 MeV, T = 270 ± 50 MeV. 



The parameters for the peak at 1.3 GeV fall between the PDG |H| values for the / 2 (1270) and 
a 2 (1320). The values for the /2(1525) are compatible with the PDG values and those for the 
structure at 2.15 GeV are compatible with both the / 2 (2150) and with the parameters for the 
structure seen in radiative J j ip decays to the same channel |14[ . 



A study has also been made of the centrally produced KgKg channel. This channel has lower 
statistics than the K + K~ channel but has the advantage that only even spins can contribute, 
which also means that there are only two ambiguous solutions to the PWA. The reaction 

pp ^ p f (K° s K° s )p s (5) 

has been isolated from the sample of events having two outgoing charged tracks plus two V°s, by 
first imposing the following cuts on the components of the missing momentum: | missing P x \ < 
20.0 GeV/c, (missing P y \ < 0.16 GeV/c and (missing P z \ < 0.12 GeV/c. 
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The quantity A, denned as A = MM 2 (pfp s ) — M 2 (KgKg), where MM 2 (pfp s ) is the missing 
mass squared of the two outgoing protons, was then calculated for each event and a cut of |A| 
< 3.0 (GeV) 2 was used to select the K$Kg channel. Requiring one V° to be compatible with 
being a Kg (0.475 < M(tt + 7t~) < 0.520 GeV) the effective mass of the other V° is shown in 
fig |5]a) where a clear Kg signal can be seen over little background. The resulting K S K S effective 
mass spectrum is shown in fig. |5|b) and consists of 2712 events. 

A Partial Wave Analysis (PWA) of the centrally produced K$Kg system has then been 
performed assuming the same analysis frame as for the K + K~ system. Since only even spins 
can contribute only S and D waves with m < 1 have been included in the PWA. The expressions 
relating the moments and the waves are given in table |2|. The PWA has been performed 
independently in 80 MeV intervals of the K° S K° S mass spectrum. 

The two complex roots, Zj, after the linking procedure, are shown in fig. ||c) and d). As 
in the case of the K + K~ channel the real parts are well separated and hence it is possible 
to identify unambiguously the two PWA solutions in the whole mass range. In addition, the 
zeros do not cross the real axis and hence there is no problem with bifurcation of the solutions. 
For one solution the mass spectrum is evenly distributed through the three D-waves and the 
S-wave is small everywhere, this solution has been ruled out. The remaining solution is shown 
in fig. ||e)-h). As in the case of the K + K~ final state, it can be seen that the S-wave shows a 
threshold enhancement; the peaks at 1.5 GeV and 1.7 GeV are interpreted as being due to the 
/o(1500) and /j(1710) with J = 0. Superimposed on the S-wave is the result of a fit using the 
parametrisation used to fit the K + K~ S-wave. As can be seen this parametrisation describes 
the K° S K° S S-wave. 

In conclusion, a partial wave analysis of the centrally produced K + K~ and KgKg systems 
has been performed. An unambiguous physical solution has been found in each channel. The 
striking feature is the observation of peaks in the S-wave corresponding to the /o(1500) and 
/j(1710) with J = 0. The D-wave shows evidence for the / 2 (1270)/a 2 (1320), the / 2 ( 1525 ) and 
the /2(2150) but there is no evidence for a statistically significant contribution in the D-wave 
in the 1.7 GeV mass region. 
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Table 1: The moments of the angular distribution expressed in terms of the partial waves for 
the K + K~ system. 



v^oo = \So\ 2 + \Po\ 2 + \Pr\ 2 + \Pi + \ 2 + \Do\ 2 + \Di\ 2 + \Dt\ 2 

V&rt w = 2\So\\Po\cos((f> s - - 4> P -) + -$-\Pq\\Dq\cos(4>p- - 4> D -) 

+ 2 -§{\P{\\D^\cos{<t> P - - 4> D -) + \Pt\\Dt\cos{<t> Pt -4> Dt )} 

v^rtn = v / 2|5' o -||Pr|co S (0 s --0 Pr )-^|Pf|| J D o |co S (0 Pr -0 D -) 

+^|p o 1l^rl^(0 P --^ r ) 
V^t 20 = ^l^l 2 -^(l^rl 2 + lA + l 2 ) + f (2|^l 2 + l^rl 2 + l^l 2 ) 

+ 2 \Sq\\ d o I cos{4> s - -4> D -) 

V^t 21 = f 5 \Pr\\Po\cos(<f> P --<f> P -) + ^\D^\\D^\cos^ D --ct> D -) 
+V2\S \\D^\cos{(/) s - - 4> D -) 

V^t 22 = $(\Pr\ 2 -\Pt\ 2 ) + $(m 2 -\Dt\*) 

= -^g{|All^r|cos(0 Pr -0 Dr ) + |P 1 + || J D 1 -|co S (0 Pi+ -0 D +)} 
+^\Po\\Do\cos(<f> P --<f> D ;) 

v^si = iglAll^o|co S (0 Pr -0 D -) + ^|P o -|Pr|co S (0 P --0 Pi -) 

V^S2 = ^{|Prl|Pr|cO S (0 p - - 4> D -) - \P+\\Dt\cOs(ct> pt - D +)} 

v^Ft 4 o = f|p»o1 2 -|(Prl 2 + Pi + l 2 ) 

>/47rt4i = ^\D Q \\D^\cos(<P D --<P D -) 
V^U 2 = >f(\Di\i-\Di\i) 
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Table 2: The moments of the angular distribution expressed in terms of the partial waves for 
the KgKg system. 



V^oo = |S "| 2 + |L>o1 2 + |L>r| 2 + |L>i1 2 
V^h = f (2|D | 2 +| J Dr| 2 + l^| 2 ) 



v^rtai = ^\Di\\Dt\cos(<f> DT - <f> D ^ 
+V2\S \\Di\cos{(/) s - -<j> D -) 

= ^(I^ri 2 -IA + I 2 ) 

y/firtu = l\Do\ 2 -l(\DT\ 2 + \Dt\ 2 ) 
V^t 41 = ^\Do\\Dt\cos(<j> D - - </> D -) 
V^t 42 = ^f(\D^-\Dt\ 2 ) 
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Figures 



Figure 1: a) The Ehrlich mass squared distribution and b) the K + K mass spectrum. The c) 
Real and d) Imaginary parts of the roots (see text) as a function of mass obtained from the 
PWA of the K + K~ system. 



Figure 2: The V^tt^lm moments from the data, are the resulting moments calculated from 
the PWA of the K + K~ final state. 



Figure 3: The physical solution from the PWA of the K + K final state. 



Figure 4: The a) S Q and b) D waves with fits described in the text. 



Figure 5: The KgKg final state, a) the M(it + 7t~) mass spectrum when the other V° is 
compatible with being a Kg. b) The K S K S mass spectrum. The c) Real and d) Imaginary 
parts of the roots (see text) as a function of mass obtained from the PWA of the K S K S system. 
e)-h) The physical solution from the PWA of the K S K S system. 
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Figure 1 
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Figure 2 
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Figure 5 
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